DYNAMIC  STRUCTURAL  MODEL 
OF  A  SUBMERGED  RING 


Jack  Thomas  Waller 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


THESIS 

DYNAMIC  STRUCTURAL  MODEL 
OF  A  SUBMERGED  RING 

by 

Jack  Thomas  Waller  Jr. 

September  1979 

Thesis 

Advisor:                  R.  E. 

Newton 

Approved  for  public  release;  distribution  unlimited 


T19196 


UNCLASSIFIED 


SECURITY   CLASSIFICATION  OF   THIS  PAGE  (Whan  Dmlm  Bnlorod) 


REPORT  DOCUMENTATION  PAGE 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.    REPORT  NUMBER 


2.  OOVT  ACCESSION   NO 


»       RECIPIENT'S  CATALOG   NUMBED 


4.      TITLE     mnd  Subtltlo) 


Dynamic  Structural  Model  of  a 
Submerged  Ring 


S.     TYPE  OF   REPORT  *  RERIOO  COVERED 

Master's  Thesis; 
September  1979 


«.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHORS 

Jack  Thomas   Waller  Jr. 


•.     CONTRACT  OR  GRANT  NUMBERS 


1      PERFORMING  ORGANIZATION  NAME  ANO  AOORE1S 

Naval  Postgraduate  School 
Monterey,  California   93940 


10.     PROGRAM  ELEMENT.  PROJECT,   TASK 
AREA  *  WORK  UNIT  NUMBERS 


11.     CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Naval  Postgraduate  School 
Monterey,  California   93940 


12.     REPORT   DATE 


September   19  79 


IS.      NUMBER  OF   PAGES 

50 


14.     MONITORING   AGENCY  NAME  *    AOORESS</f  Itttmrmnt  Iron   Controlling  Otlico) 


IS.     SECURITY  CLASS,  (ol  thi,  riport) 

Unclassified 


t»«.     DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


16.     DISTRIBUTION   STATEMENT  (ol  thlm  Rmpmtt) 


Approved  for  public  release;  distribution  unlimited 


17.     DISTRIBUTION  STATEMENT  (ol  tho  ftroct   mirnd  In  Btoek  30,  II  altlotont  from  Rtpoti) 


IS.     SUPPLEMENTARY  NOTES 


It.     KEY  WORDS  (Conrimi.  on  rmrota 


tarr  mnd  Identity  or  bloc*  nimnmor) 


Structural  Ring 
Trigonometric  Series 
Underwater  Shock 
Fluid  Structure  Interface 


20.     ABSTRACT     Conttnuo  an  <■«»•»■« 


•-mt  mnd  Idmniltr  or  mlomk 


toot) 


A  dynamic  structural  model  of  a  submerged  ring  is 
developed  using  trigonometric  series.  It  is  constructed 
for  use _ in  conjunction  with  a  finite  element  fluid  model 
to  examine  the  effects  of  cavitation  on  underwater  shock 
loading  of  a  structure.  The  governing  equations  and  the 
time  integration  algorithm  used  in  the  model  are  described 


DO  ,  :°:M7J  1473 
(Page  1) 


EDITION  OF    I  NOV  St  IS  OBSOLETE 
S/N    0102-014-6601   I 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAOE  (Whon  Dmlm  Bnlmrmd) 


UNCLASSIFIED 


t»euwTv  ci*m»ic*TiQM  qw  run  'tat^xi  n»<«  t« 


(20.   ABSTRACT   Continued) 


Results  predicted  by  the  model  are  compared  to  known 
results.   The  program  listing  is  given. 


DD     Form       1473 

,  1  Jan  73  ,,  UNCLASSIFIED 

S/N    0102-014-6601  2 


||CU*lty  eiAMI'lCATlON  0'  THI$  **0*C»*««  0«««  *««•»•*) 


Approved  for  public  release;  distribution  unlimited, 

Dynamic  Structural  Model 
of  a  Submerged  Ring 

by 

Jack  Thomas  jailer  Jr. 
B.S.M.E.,  New  Mexico  State  University,  1975 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  MECHANICAL  ENGINEERING 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
September  1979 


NAVAL  POSTGRADUATE  SCHOOl 
MONTEREY,  CA  93940 


ABSTRACT 

A  dynamic  structural  model  of  a  submerged  ring  is 
developed  using  trigonometric  series.   It  is  constructed 
for  use  in  conjunction  with  a  finite  element  fluid  model 
to  examine  the  effects  of  cavitation  on  underwater  shock 
loading  of  a  structure.   The  governing  equations  and  the 
time  integration  algorithm  used  in  the  model  are  described 
Results  predicted  by  the  model  are  compared  to  known 
results.   The  program  listing  is  given. 
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I.   INTRODUCTION 

The  program  evolved  during  this  study  models  a  submerged 
circular  ring  using  trigonometric  series.   It  was  developed 
for  use  in  conjunction  with  a  finite  element  fluid  model 
based  on  a  displacement  potential  formulation.   The  purpose 
of  the  combined  models  is  to  predict  the  effects  of 
cavitation  on  underwater  shock  loading  of  the  structure. 
More  information  on  the  fluid  model  can  be  found  in 
References  1  and  2.   The  purpose  of  this  paper  is  to 
describe  the  development  of  the  structural  model  and  to 
present  some  of  the  results  obtained  from  its  use  in 
combination  with  the  fluid  model.   A  listing  of  the  FORTRAN  IV 
program  implementing  the  structural  model  is  given  in  the 
Appendix. 
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II.   GOVERNING  EQUATIONS 

A.   DIFFERENTIAL  EQUATIONS  FOR  THE  DEFLECTION 
OF  A  THIN  CIRCULAR  RING 

In  Figure  1  the  solid  line  represents  the  ring  after 

deformation,  and  the  dotted  line  the  ring  before  deformation 

For  small  deflections  the  curvature  of  the  element  nun.  can 

be  taken  as 


_1_  =   d9  +  Ad9 
R,     ds  +  Ads 


where  R,  w,  9  and  s  are  defined  as  shown  in  Figure  1  and 
w  is  taken  positive  inward.   Making  use  of  the  relations 


2  2 

ajo      jn     jn      dw    dw,  dw      dw. 

Ad9   =   d9,  -  d9   =   3—  +  — *-  ds  -  -g—  =   — x-  ds 

1          ds    ,  2  ds      ,2 

ds  ds 

ds 
Ads   =   ds   -  ds   =   (r-w)de  -  rd9   =   -wd9   =   -w  — 


and  substituting  into  the  equation  above  yields 

2  2 

d9  +  *«   ds  ,.n  .  w.      14  ds(l  +  |) 

i   =        ds2  d6(1  +  R}      ds2         R 

K_        —    T 


ds(l-|)        ds(L-|)(l  +  g)    ds(l-|)(l  +  |) 


Neglecting  the  higher  order  terms,  the  above  expression 
becomes 
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d(jj~  ,  dor  j- 


Figure  1.   Geometry  of  Ring  Deformation 
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2  2 

57  "     di(1  +  R)+— 2     -     r(1  +  r)+TT 

1  ds  ds 


or   alternately 


2  2 

11      w,dw      1,  dwx  ,, » 

R7  "  R   -  12  +  TT  =  -2  (w  +  — 2>             (1) 

1          R     ds       R  d0 


For  a  ring  where  the  thickness  is  small  compared  to  the 
radius  and  elastic  behavior  is  assumed,  it  can  be  shown 
that  the  approximate  relationship  between  deflection  and 
loading  is  [Reference  7] 


i  -  I  -   -  *  (2) 

R„    R        D  K    ' 


where  M  is  the  bending  moment  about  the  centroidal  axis 
and  D  is  the  flexural  rigidity  of  the  ring.   A  positive 
bending  moment  produces  compression  in  the  outside  fibers 
of  the  ring.   Combining  equations  1  and  2  yields  the 
differential  equations  for  the  deflection  of  the  ring  given 
below. 


—7T    ( j  +   w)   =   -  -  (3) 

R    99^  U 
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B.   GOVERNING  EQUATIONS  OF  MOTION 

The  governing  equations  of  motion  used  in  the  program 
were  arrived  at  by  the  application  of  Hamilton's  principle. 
In  order  to  apply  Hamilton's  principle,  it  is  first  necessary 
to  derive  expressions  for  the  strain  energy  and  kinetic 
energy  of  the  ring  as  well  as  the  work  done  by  the  external 
loads.   In  these  derivations  the  ring  was  taken  to  be 
homogeneous,  elastic,  and  of  unit  width.   The  pressures  on 
the  ring  and  its  deflection  were  represented  by  the 
pressures  at  and  the  deflection  of  a  set  of  nodal  points 
equally  spaced  along  the  circumference  of  the  ring. 

.  The  shock  front  is  assumed  to  approach  the  ring  normal 
to  the  6=0  plane  as  shown  in  Figure  2,  where  9  is  taken 
to  be  positive  counterclockwise. 


Figure  2.   Ring,  Shock  Front  Orientation 
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Since  the  ring  is  symmetric  about  the  8=0  plane, 
results  are  given  for  9  from  0°  to  180°.   The  response  of 
the  ring  is  represented  by  the  response  of  a  set  of  N+l 
nodal  points  equally  spaced  along  the  ring.   The  radial 
displacement  is  approximated  by  a  trigonometric  cosine 
series 


N 
w  =    I        a   cos  n9    .  (4) 

n=0    n 


The  tangential  deflection  is  represented  by  the  trigonometric 
sine  series 

N 

v  =    y   b   sin  n9    ,  (5) 

K    n 
n=l 

where  v  is  taken  to  be  positive  in  the  negative  9  direction. 
Similarly  the  normal  pressure  applied  to  the  ring  is 
represented  by 


N 

J    c   cos  n9   .  (6) 

n=0 


1.   Strain  Energy 

The  strain  energy  is  comprised  of  two  components 
[Ref .  4] .   One  component  is  the  strain  energy  due  to  bending 
and  the  other  to  strain  in  the  9  direction  (£Q)-   The  strain 
energy  stored  in  the  ring  due  to  bending  is  calculated  as 
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the  work  done  by  the  internal  moment  acting  over  the 
entire  ring,  and  for  a  ring  of  unit  width  is 


2tt 
UB   =   f     *       M  Rlcd6 


where  k  is  the  curvature  of  the  ring.   Combining  this 
equation  with  equation  3  and  using  the  relation  <  =  M/D 
yields 

,     2tt    .   .2      ~ 

uR  =  i   /   d  i  (1-g  +  w)2  de  (?) 

o    R   3e 

The  strain  energy  due  to  the  strain  in  the  9  direction  is 
obtained  as  the  distance  traveled  by  the  average  normal 
force  in  the  9  direction. 

The  circumferential  strain  is  composed  of  two 
components  the  first  of  which  is  the  result  of  the  radial 
displacement  of  the  ring  ( -5-  )  and  the  second  results  from 
tangential  displacement  and  is  (—  -*-r)  .      The  strain  energy 

R  do 

from  the  normal  force  (N)  is 

,  2tt  ,  2tt 

UQ      =      ^      /         NeeR   d9      =      j      j         A£QE£eR  d9 

0  0 

or  (8) 

1  2tt  _  ,  2tt    ,  «  - 

U,       =      y      /         RAE£.Z    d9      =      4      /         ^f    (|^  +    w)Z    Rd9 

2  0  2       0         R2       39 
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where  A  is  the  cross  sectional  area  of  the  ring  and 

E  is  the  modulus  of  elasticity  of  the  ring  material. 

Combining  equations  7  and  8,  the  expression  for  the 
total  strain  energy  is 

1       f2lT     D      ,82w  ,2     ,.  1       r2Tr   AE,8v  .2     ,_ 

ut  -  2  /   71  (^2  +  w)  de  +  I  /   -R-(ae  +  w)  de 

o   R   3e  0 

u  u  (9) 

Substituting  into  this  equation  the  series  expressions  for 
radial  and  tangential  displacements,  and  assuming  constant 
geometrical  and  material  properties  over  the  ring  yields 
after  integration 


N  N                * 

TTD  T      "  2,2   .,2,  AEtt  .5   ,  .        ,2, 

UT  =   773  [  L  *n  (n  "  1)  +  "2R-  I  I       (nbn  +  a  )  ] 

2R     n=0  n=0 


An  expression  for  the  change  in  strain  energy  5U 
corresponding  to  small  changes  in  tangential  and  radial 
displacements  (6b  and5a  )  can  now  be  found 

su  =  ID  f  f   (n2.1)2       +  MI  [  f   (n\  +  nan)5bn 
n=0  n=0 

+  (nb   +  a)Sa)]       (10) 
n    n    n 


* 
In  this  equation  and  several  that  follow  which  contain 
a  summation  expression  starting  with  n  =  0,  a  factor  of  2 
associated  with  the  n  =  0  term  is  omitted  for  the  sake  of 
brevity.   The  term  is  accounted  for  in  the  solution  process 
used  with  the  model. 
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2.   Work 

The  work  associated  with  a  small  change  in  radial 
displacement  is  given  by 


2ir 
6W  =    /    p  6w  R  d6 
0 


Substituting  the  relations 


and 


N 
p  =      c   cos  m( 
m=0 


N 

5w  =    7   5a   cos  nf 

n     n 

n=0 


into  this  expression  yields,  after  integration 


N 


<SW  =    J   c   5a   ttR 

n       n    n 

n=0 


3 .   Kinetic  Energy 

The  expression  for  the  kinetic  energy  of  the  ring 
can  be  arrived  at  by  considering  an  infinitesimal  section 
of  the  ring  that  has  been  set  in  motion.   The  kinetic  energy 
of  the  element  can  be  represented  as  the  sum  of  a  contribu- 
tion due  to  translation  of  the  mass  center  and  a  contribution 
due  to  rotation.   This  leads  to 
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dT   =   |  p(ARd9)  (w2  +  v2)  +  |  pIRd6(^-)2 


where  I  is  the  area  moment  of  inertia   about  the  centroidal 
axis. 

Using  the  series  relations 


N  N 

w  =    7    a   cos  n6        w  =    Y    a   cos  n6 
n=0  n=0 


N  N 

v  =  b   sin  n0        v  =    Y   b   sin  n( 

_-  n                      -i    n 

n=  1  n=l 


yields 


1  N   •         ?     N  9 

dT  =  j   p (ARd9 ) [ (  I      an  cos  n0 )   +  (  £   b   sin  nG)  ] 

n=0  n=l 


1  1   N  2 

+  i  pIRd6(-  ~  I      a_  n  sin  n6) 

^  r\     ■»    xi 

n=l 


Integration  over  the  ring  yields 


T   =   *£       [ (AR  +  n2  h a  2  +  ARb  2 
2   n=0  R   n 


The  resultant  change  in  kinetic  energy  due  to  a 
small  change  in  velocity  components  can  then  be  found  as 


N  -  .   . 

6T   =   ttp   J   [  (AR  +  n   £)a  6a   +  ARb  5b  ]         (12) 
L  rt  R  n   n      n   n 

n=0 
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4 .   Application  of  Hamilton's  Principle  To  Find 
Coupled  Equations  of  Motion 


From  Hamilton's  principle  [Ref.  5]  it  is  known 


that 


/ 

t. 


(5T  -  5U  +  5W)  dt   =   0 


(13) 


Combining  equations  10,  11,  12  and  13  yields,  after 
carrying  out  an  integration  by  parts  on  the  first  term, 


N 


2  I, 


I        ttP(AR  +  n   4)  ^r,6^ 
n=0  R   n   n 


2   .. 

f    a   6a   dt] 
{.     n   n 


N 
+        I        TrpAR[b      5b 
nio  n        n 


/  b      6b      dt] 

I  n        n 

rl 


VttD    r    ?    ,2     ,,2 


N 


/  '{™   [     £     (n2-l)"    a      5a    ]    +  £f!   [      J       (n2b      +   na    )  6b 

L3"  nn  ^  n  n 


ti     R         n=0 


n        rr 


n=0 


n        n 


N 

+    (nb      +   a    ) 5a    ]    -      T        c      5a      ttR}        dt      =      0 
n  n        n  Ln        n        n 

n=0 


Since   the    5a      and    6b      can   be   chosen   arbitrarily/    this 

n  n  J 


becomes 


,,  ,r,2      2    **  ,c>2r,  ,rx2.    2      ,.2,  ,c>  2  n 

(1    +    (R}       n    )an   +     (R}     [1    +     (R>     (n     "1}     ]an   +     (R}         =      pA 


*n  +    (I)2    nan+    (l)2    n\      -      ° 
1/2 


(14) 


where   c  =    (E/p)  is   the    speed  of   sound   in    the   ring   and 
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where  r  is  the  centroidal  radius  of  gyration  of  the  ring 

2 
section  and  the  relation  I  =  Ar   has  been  used. 

5 .   Separated  Equations  of  Motion 

It  is  common  practice  to  separate  ring  deformations 

into  flexural  modes  and  extensional  modes.   The  initial 

version  of  this  program  utilized  such  a  resolution.   The 

flexural  mode  is  characterized  by  the  requirement  that 


eQ  =   "  (w  +  g^)/R  =   0  . 


This  implies  that  the  Fourier  coefficients  satisfy  the 
relation 


a     +  nb      =0  (15) 

n(F)      n(F) 


The  extensional  mode  is  defined  to  be  geometrically  orthogonal 
to  the  flexural  mode  so  that 


na      -  b      =0  (16) 

n(E)     n(E) 


Substituting  equation  15  into  equations  14  yields  the  result 


2 
[pA((n  t1)  +  zn2)]^     +  [^z(n2-l)2]  a^     =   cn 
n2  n(F)     R2  n(F)       n 

b^     =  -a    /n  (17) 

nCF)      nCF) 

2 
where  z  =  (r/R) 
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If  equation  16  is  substituted  into  equations  14,  the  result 
is 


[pA((n2  +  l)  +  zn2)]a     +  [^  (  (n2  +  1)  2  +  z(n2-l)2]a 

n(E)     R  n(E) 

"   cn  (18) 


b      =   na 
n(E)        n(E) 


Initial  solutions  were  carried  out  using  equations 
17  and  18  and  then  combining  the  results  by 


a   =   a     +  a 

n(F)     n(E) 


b   =   b     +  b 
n       n(F)     n(E) 


This  procedure  gave  results  which  satisfactorily  predicted 
the  ring  bending  moments  but  gave  poor  accuracy  for  axial 
force  determinations. 

Re-examination  of  the  foregoing  solution  process 
revealed  that  the  mode  shapes  defined  by  equations  15  and 
16  are  not  orthogonal  with  respect  to  the  mass  or  the 
stiffness  matrices  of  the  system.   It  was  accordingly 
decided  to  return  to  equations  14  and  solve  them 
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simultaneously  in  the  time  integration  algorithm.   This 
change  led  to  accurate  axial  force  determinations .   The 
separated  equations  of  motion,  17  and  18,  are  used  in  the 
program  as  described  in  Section  III,  b. 
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III.   TIME  INTEGRATION  METHOD 

A.   CENTRAL  DIFFERENCE  ALGORITHM 

Time  integration  is  accomplished  using  a  central 

difference  algorithm.   If  a  ,  b  ,  and  c   are  known  at 
^  n   n       n 

time  t(   ,  then  a(l)  (the  value  of  a   at  t(l))  may  be 

n  n        i        j. 

evaluated  from  the  first  of  equations  14.   Letting  h 

.  (i-k) 
represent  the  time  step  and  a       represent  a   at  time 

(i) 
t    -  h/2,  the  next  value  of  a   is  calculated  from 


kU+H)     m     •  (i-H)   +  h»(i)  (19) 

n  n  n 


Using  this  result,  the  next  a   (at  time  t      =  t    +  h) 

n 

is  calculated  from 


a(i+1>   =   a(i)  +  h  a(i+l5»  (20) 

n  n        n 


The  value  b     is  similarly  found  from  the  second  of 

n  1 

equations  14.   A  pair  of  equations  paralleling  19  and  20 

.  ( i+i^)       h  ) 
are  used  to  calculate  b    2   and  b    .   When  these  steps 

n  n  r 

have  been  completed  for  each  n,  the  value  of  radial 
displacement  w  is  found  at  each  structural  node  from 
equation  4.   These  displacements  are  passed  to  the  fluid 
program  where  they  furnish  required  interface  boundary 
conditions  to  allow  an  advance  to  time  t      .   Values  of 
fluid  pressure  at  the  interface  nodes  are  returned  by  the 
fluid  program. 
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The  new  nodal  pressures  are  used  to  calculate  c  ' s  at 

r  n 

t      .   It  is  then  again  possible  to  advance  a„  and  a 

n      n 

using  equations  14,  19,  and  20  and  to  perform  the  parallel 

calculations  for  b   and  b  . 

n      n 

The  solution  process  is  started  with  the  ring  at  rest 

under  uniform  hydrostatic  pressure.   Under  the  loading  a- 

•  (0) 
is  the  only  nonzero  Fourier  coefficient.   Because  a    =  0 

1  n 

is  known  (rather  than  a     ,  a  fictitious  starting  value  of 

a     is  first  calculated  from 
n 


•(-%)   _    h  »(0) 

n  2   n 


A  similar  starting  procedure  is  used  for  b  . 

B.   SELECTION  OF  THE  NUMBER  OF 
VIBRATIONAL  MODES  USED 

The  number  of  vibrational  modes  that  can  be  modeled 
accurately  is  limited  by  the  numerical  integration  algorithm. 
Specifically,  the  algorithm  becomes  unstable  for  time  steps 
in  excess  of  about  0.3  of  the  period  of  the  structural  mode. 
The  accuracy  of  the  algorithm  deteriorates  even  before  the 
stability  limit  is  reached.   For  this  reason,  a  criterion 
was  established  for  limiting  the  number  of  modes  used, 
based  on  the  time  steps  selected. 

In  the  program  the  separated  equations  of  motion  for 
the  extensional  and  flexural  modes  were  used  to  find  the 
frequencies  needed  in  applying  the  criterion.   From 
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equations  17  and  18  the  natural  frequencies  of  the 
extensional  and  flexural  modes  were  found  to  be 


u>  2     =   (§)2(1  +  n2)  n  =  0,1,2,3,... 

n(E)       R 

2   2     2 
2        -c,  2  ,r,  2  n  (n  -  1)  -^    ->    A 

u>      =   («)  (p)   S —    n  =  2,3,4,.. 

n(F)       R    R     n2  +  l 


The  criterion  used  was  that  the  period  for  the  highest 
mode  retained  be  at  least  five  times  the  time  step  used. 
Since  the  extensional  modes  have  higher  frequencies  for  a 
given  n,  this  criterion  comes  into  effect  first  for  the 
extensional  modes.   Once  the  point  is  reached  (where 
n  =  n_ (max) ) ,  the  program  switches  from  the  coupled  equations 
of  motion  14,  to  the  separated  equation  17.   When  this 

occurs  the  higher  extensional  mode  coefficients  (n  >  n_(max)) 

hi 

are  not  needed,  and  only  the  flexural  mode  coefficients 

(for  n  >  n_ (max) )  are  used.   If  the  time  step  is  large 

enough  that  the  criterion  is  also  met  by  the  flexural 

modes  (n  =  n_ (max) )  then  both  the  a  's  and  b  's  are  omitted 
F  n        n 

for  all  modes  where  n  is  greater  than  n_ (max) . 

r 
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IV.   RESULTS 

To  check  the  accuracy  of  the  results  obtained  from  the 
program,  two  kinds  of  comparisons  were  made.   First,  the 
solution  of  equations  14  arrived  at  by  the  program  were 
compared  to  an  analytic  solution  of  equations  14.   Second, 
to  verify  that  equations  14  correctly  predict  ring  behavior, 
calculated  dynamic  response  of  the  ring  to  two  particular 
loadings  was  compared  to  known  static  results  [Reference  6] 
for  the  same  loadings . 

A.   COMPARISON  OF  MODEL  RESULTS  TO  EXACT  SOLUTION 

The  check  on  solution  accuracy  was  based  on  the  exact 
solution  for  a  ring  initially  at  rest  and  suddenly   loaded 
by  a  steady  pressure 

P   =   [1.625  -  .25  cos  0  -  .375  cos  26]  MPa 

The  exact  solution,  using  the  parameters  specified  in 
Table  I,  is 


aQ   =   4.05844xl0"3 [1  -  cos  (1011. 303t) ] 

a,   =   -159.483t2  -  1.5617xl0"4  [1  -  cos  (1429. t) ] 

a2   =   -.104129[1  -  cos  (85.616)t] 

-  3.7355xl0~5[l  -  cos  (2259. 7t)]        (22) 


Table  I  gives  results  for  radial  displacements  at 
t  =  1  msec  for  0  =  0°,  90°  and  180°.   Computer  results  are 
given  for  four  different  time  steps.   Even  the  coarsest 
time  step  gives  results  within   1%  of  the  exact  values 
obtained  from  equations  22. 

B.   COMPARISON  OF  MODEL  RESULTS  TO  STATIC  BEHAVIOR 

To  compare  the  dynamic  response  predicted  by  the  program 
to  the  known  static  response  of  the  ring  for  the  same 
loading,  two  simple  cases  of  loading  were  assumed.   The 
first  loading  case  examined  was  uniform  pressure  surrounding 
the  ring;  the  second  case  was  two  equal  and  opposite  concen- 
trated loads  acting  180°  apart  as  shown  in  Figure  3.   The 
second  loading  case  could  not  be  represented  exactly  by 
the  finite  trigonometric  series  loading  representation  used 
by  the  program.   It  was  approximated  by  nonzero  values  of 

pressure   only  at  8  =  0°  and  8  =  180°. 


B 


Figure  3 .   Concentrated  Loading  on  Ring 
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TABLE    I 


Time 
Step 

sec 

AR 
NPl 

mx  10~7 

AR 
NP2 

mx  10"7 

AR 
NP3 

mx  10"7 

1 

AR 

predicted 
by  program 

2.5xl0"4 

11619. 

23687. 

17556. 

LxlO"5 

11679 

23472. 

17553 

LxlO'6 

11682 

23469. 

17555. 

LxlO"7 

11683. 

23469. 

17555. 

AR  from 
equations  17 

11683. 

23469. 

17555. 

Ring  Parameters: 


R  =  5  m,  A  =   .05  m  ,  E  =    200  GPa 
p  =  7830  kg/m3,  r  =   .158m 


Comparison  of  program   solution  with 
analytic   solution 


TABLE    II 


AR    (max) 


mm 


Axial 
Force 
(Max) 

MN 


Program  results 


4.98 


-9.8 


Twice  Known 
Static  Results 


5.00 


-10.0 


Ring  Parameters: 


R  =  5  m,  A  =  .05  m  ,  E  =  200  GPa 
p  =  7830  kg/m3,  r  =  .158  m 


Comparison  of  program  dynamic  solution 
with  twice  known  static  solution  for 
uniform  external  pressure  on  ring 
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Comparison  between  the  dynamic  and  static  responses 
of  the  ring  were  made  for  the  radial  displacements,  bending 
moments/  and  forces  at  particular  nodal  points  on  the  ring. 
For  the  uniform  pressure  case  any  nodal  point  on  the  ring 
can  be  used  since  the  displacements  are  uniform  over  the 
ring.   For  the  concentrated  loading  case,  the  nodal  points 
chosen  were  at  the  points  of  application  of  the  forces  and 
midway  between  the  application  points. 

For  the  dynamic  response  of  the  ring  predicted  by  the 
program,  the  ring  is  taken  to  be  initially  at  rest  and 
undeformed.   The  anticipated  behavior  of  the  ring  after 
the  application  of  the  loading  is  for  the  ring  to  deflect 
through  the  radial  displacements  found  for  the  static 
loading  and  to  continue  to  deflect  until  a  radial  displace- 
ment with  a  magnitude  of  approximately  twice  the  static 
displacement  is  reached.   The  maximum  deflections  were 
expected  to  occur  after  a  time  approximately  equal  to 
one  half  the  period  of  the  vibrational  mode  that  dominates 
the  deflection  of  the  ring.   For  the  uniform  loading  case 
the  dominant  vibrational  mode  is  the  fundamental  extensional 
mode.   For  the  concentrated  loading  case,  it  is  the  second 
flexural  mode.   The  periods  for  these  vibrational  modes  can 
be  found  using  equations  21.   Tables  II  and  III  give  the 
comparisons  between  the  dynamic  response  from  the  program 
and  twice  the  known  static  response  for  the  two  loading 
cases  considered.   As  is  seen,  relatively  good  correlations 
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TABLE  III 


Slumber  of 
SJodal  Points 
Used 

i —  —  — 
AR^max) 

m 

AR^(max) 
m 

Axial 
Force 
at  B 

MN 

Bending 
Moment 
at  B 

MN»m 

5 

Model 
Results 

.2965 

-.2839 

-3.760 

-7.910 

Static 
Analysis 
Results 
Times  2 

.2634 

-.2421 

-3.536 

-6.434 

7 

Model 
Results 

.2013 

-.1857 

-2.353 

-4.737 

Static 
Analysis 
Results 
Times  2 

.1862 

-.1713 

-2.500 

-4.543 

9 

Model 
Results 

.1510 

-.1393 

-1.743 

-3.544 

Static 
Analysis 
Results 
Times  2 

.1425 

-.1311 

-1.913 

-3.477 

Ring  Parameters: 


R  =  5  m,  A  =  .05  m  ,  E  =  200  GPa 
p  =  7830  kg/m3,  r  =  .158  m 


Comparison  of  dynamic  solution  with  twice  known 
static  solution  for  loading  of  Fig.  3 


32 


between  the  dynamic  responses  and  twice  the  static  responses 
were  obtained.   The  maximum  values  shown  in  Table  III  occur 
at  approximately  35.3  milliseconds  after  time  zero,  which 
compares  favorably  to  one  half  of  the  period  of  the  second 
flexural  mode  (36.6  milliseconds).   The  maximum  values  in 
Table  II  occur  at  3  milliseconds  which  is  approximately 
one  half  the  period  of  the  fundamental  extensional  mode 
(3.11  milliseconds). 

C.   OUTPUT  EXAMPLES 

A  sample  of  the  output  produced  by  the  structural 
program  is  shown  in  Figures  4  and  5.   The  print  code  at  the 
top  of  Figure  4  is  an  input  code  used  to  select  the  informa- 
tion desired  in  the  print  out.   Following  the  print  code 
are  two  sets  of  input  parameters  to  be  used  in  the  structural 
model.   These  are  followed  by  the  mass  and  stiffness 
coefficients  calculated  for  each  of  the  vibrational  modes . 
This  information  remains  constant  for  any  given  run  and 
is  therefore  only  printed  once.   The  rest  of  the  information 
varies  with  time  and  can  be  printed  as  often  as  desired. 
The  pressures  at  nodal  points  generated  by  the  fluid  model 
are  printed  out;  these  pressures  are  for  one  time  step 
before  the  time  printed  below  the  pressures.   Following  the 
time  come  the  Fourier  coefficients  for  the  radial  displace- 
ments, tangential  displacements,  and  pressures.   In 
Figure  5  at  the  top,  the  bending  moment  and  the  axial  force 
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at  each  nodal  point  are  printed  for  the  time  given  in 
Figure  4.   Tabulation  of  the  radial  displacement  follows. 

Figure  6  is  an  example  of  the  graphical  output  obtained 
from  the  fluid  model  when  the  fluid  and  structural  models 
are  run  jointly.   It  is  a  time  sequence  of  eight  plots  of 
fluid  nodal  pressures  over  the  domain.   The  plots  are  shown 
at  8  ms  intervals.   The  left  hand  edge  of  each  plot  repre- 
sents the  plane  of  symmetry  and  the  top  and  bottom  rows 
represent,  respectively,  the  entry  and  exit  faces  for  the 
shock.   The  x's  on  the  lower  left  side  of  the  plots 
represent  dummy  nodes  inside  the  structure.   The  pressure 
ranges  for  the  mapping  characters  are  shown  in  Figure  7 . 
In  Figure  6  the  development  of  the  cavity  can  be  followed. 
It  develops  in  the  second,  third,  and  fourth  frames, 
collapses  in  the  fifth,  and  has  vanished  by  the  sixth. 
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PRINT  CODE:  121111111  1  lie  11415161718  1  1 
NNEL  NNP   M    N  NNPftlN  NNPflSM 


16 


17    2    @    '?   17 


fl        E        R        RI 
500D-9 1  @ .  20@E  12  3 .  66     0 . 1 58 


KKflSN  XMfilM 

0.0  391.50 

0 .  1 02 1 7Ii  07  39 1 .  39 

0.20434D  07  393.07 

y.3065 II'  07  395.02 

0.40S6SD  07  397.76 

0.510S6D  07  401.29 

0.61303D  07  405.59 

0.71520D  07  410.63 

0.317371!  07  4-16.56 

0.91954D  07  ^23.21 

0.102171!  08  420.65 

0.11239D  38  438.87 

0.1 226 ID  08  447.38 

0. 132S2D  08  457.66 

0 . 1 4304D  08  463 . 23 

0.15326D  03  478.59 

0. 16347D  OS  491.72 


I 

XKfllN 

1 

0. 

40040D 

09 

2 

0. 

40  000  D 

09 

3 

0. 

403t0D 

09 

4 

0. 

42560D 

09 

5 

0. 

49000D 

09 

6 

0. 

63040D 

09 

7 

0. 

39000D 

09 

8 

0. 

1321 6D 

10 

9 

0. 

19876D 

10 

10 

0. 

29660D 

10 

11 

0. 

43204D 

10 

12 

0. 

61600D 

10 

13 

0. 

35796D 

10 

14 

0. 

11690D 

11 

15 

0. 

15610D 

11 

16 

0. 

80470D 

11 

17 

0. 

2641 0D 

11 

RHO 

0.733D  04 
XKE1N 

0.0 

0.40000D  09 

0 . 30000D  09 

0. 12000D  10 

0. 16000D  10 

0.20000D  10 

0 . 24000D  1 0 

0.23000D  10 

0.32000D  10 

0.36000D  10 

0.40000D  10 

0. 4400SD  10 

0.48000D  10 

0.52000D  10 

0.56 000 D  10 

0.60 000 D  10 

0.64000D  10 


NODRL  POINT  PRESSURES 


0. 12D  08  0.98D  07  0.56D  07  0.2  ID  07 
0.21D  06  0.2OD  06  O.20D  06  0.2OD  06 
0. 16D  06  0.20D  06  O.20D  06  O.20D  06 


0.52D  06  0.23D  36 
0.21D  06  0.21D  06 
0.2OD  06 


TIME  =   0.250  MS 


EN 

0.  0 

■0. 43636D-05 
■0. 76030D—05 
-0.9095211-05 

-0.37862D-05 

-0.70528D-05 

-0.48681D-05 

■0.27812D-05 

•0. 11520D-05 

-0. 24765Ii-06 

0.76762D-O3 

0.70517D-07 

O.  11450D-07 

O.64630D-07 

0.48831D-07 

0.38325D-07 

0.33331D-07 


AN 

O. 59923 D-03 
0.19121D-03 
0. 16725D-03 
3. l34c5B-03 

q  .  £  39 1 S  B— 04 

0. 69591 B-05 
0.22223D-05 
-0. 76762D-07 
-0 .  77569D— 06 
-0. 13740B.-06 
-0.8461  SB— 06 
-0.68433B-06 
■0 .  5748715—06 
•0.54130D— 06 


CM 

3.  16297TJ  07 
0.27357B  07 
0.23797D  07 

IS.  Ico'^lD  0r' 
O. 13588D  07 
6.864201!  36 
0.48802D  06 
0.231 7211  06 

76529. 

18334. 
-8451.9 
-15526. 
--825.4 
-14242. 
-10768. 
-8187. 1 
-*1  :-;•=•  7 


Fig.  4.   Output  Exacple 


35 


NP 


MOMENT 


FORCE 


6 


1@ 

11 

12 
13 
14 
15 
16 
17 


61466. 

I  -I ;-'  i-  i^. .-' 

-46163. 
-37393. 
-7463. 4 
-5081.3 
-5665. 4 
-5152.0 
-4923.6 

-3399. 1 
-7444.2 
-3635.2 
-5265. S 
-4963.3 
-4983. 9 


y.-~ 

.227171: 

07 

V  ' 

.21430D 

07 

0 

1  c  i  Cr~'Tj 

i"j~ ? 

0 

!l4334D 

07 

0 

. 1144SB 

07 

@ 

. 10173D 

07 

0 

. 9 9 940 B 

06 

0 

. 1 0045D 

07 

u 

, 10S06D 

07 

Pi 

.  9982513 

06 

i" 

.  100171! 

07 

0 

. lOOlOD 

07 

0 

.  yy?yc!D 

Ob 

0. 

. 99765D 

06 

0  ■ 

. 10009D 

07 

0i 

.  99989D 

06 

0, 

.  99806D 

06 

NP 


MODAL  POINT  DISPLACEMENTS 


p; 


1 

0. 13157H-02 

c 

0. 11734D-02 

3 

0.36677II-03 

4 

0 . 62037D-03 

5 

0.51322D-03 

6 

0 . 50099D— 03 

r 

0.50035D-03 

.". 

0.49977D-03 

q 

0. 4998 ID— 03 

10 

O . 49984D-03 

11 

0 . 50032D-03 

12 

0. 50019D-03 

13 

0. 49695D— 03 

14 

0. 49939D— 03 

15 

0.49947D-03 

16 

0.49921D-03 

17 

0 . 49944D-03 

T=10. 

40---12.73 

10.23.52 

Fig.  5.   Output  Example  (cent.) 
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Symbol  Pressure  Range  (m  pa) 
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Figure  7.   Pressure  Code 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  comparisons  made  between  the  program  results  where 
convergence  of  the  solution  has  been  obtained  to  results 
obtained  by  other  methods,  show  that  the  model  satisfactorily 
predicts  the  dynamic  response  of  a  ring  to  external  pressures 
The  ring  model  cannot  predict  exactly  the  behavior  of  a  real 
structure  such  as  a  submarine  hull.   However,  by  a  suitable 
choice  of  ring  parameters  the  response  of  the  ring  structure 
will  provide  a  good  first  approximation  to  the  response  of 
a  real  structure. 

A  useful  improvement  to  the  program  would  be  to  include 
a  section  which  would  automatically  pick  out  maximum  values 
of  bending  moments  and  axial  forces  in  the  ring  as  well  as 
when  and  where  they  occur.   This  would  allow  any  desired 
time  intervals  between  data  printout   to  be  used  without 
missing  peak  responses. 
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